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Abstract. In this paper we classify the free boundary associated to equihb- 
rium configurations of compressible, self-gravitating fluid masses, rotating with 
constant angular velocity. The equilibrium configurations are all critical points 
of an associated functional and not necessarily minimizers. Our methods also 
apply to alternative models in the literature where the angular momentum per 
unit mass is prescribed. The typical physical model our results apply to is that 
of uniformly rotating white dwarf stars. 

1. Introduction 

In this paper we study the free boundary associated to rotating star models of 
white dwarf stars with prescribed constant angular velocity. Thus we are consider- 
ing figures of equilibrium for compressible, self-gravitating fluid masses. 

There has been a tremendous amount of work on incompressible, self-gravitating 
fluid masses rotating with prescribed constant angular velocity since the primary in- 
vestigations by Newton. Various mathematicians like MacLaurin, Jacobi, Dirichlet, 
Riemann, Poincare, H. Cartan and Chandrasekhar made significant contributions 
to the field, studying bifurcation sequences and analyzing the stability of various 
equilibrium shapes. A historical account and details of these investigations may be 
found in Chandrasekhar 's treatise [10] and Tassoul's book [17]. 

In the compressible case for the model with prescribed constant angular velocity 
w > (cf. [14]), we consider the functional 

(1.1) J{p)= f A{p)d^-\( io'r^pd^- [ pBpdi, 

where ^ = (6,6, 6), r = ^M+^, P>0,P^ i'(R') n L°° (M^) and 
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Moreover we impose the constraint 

(1.2) I P=l. 

As p(^) represents the density of the stellar material, (1.2) means that the mass 
of the star is prescribed. Later in our paper we will assume in addition that the 
density p(^) is axisymmetric, i.e. 

P(0 = p{r, z) where z{^) = ^3. 

This assumption means that the star is rotating about the ^3-axis. The function 
A{p) is the pressure and thus represents the equation of state of the stellar mate- 
rial. We assume that A{p) E C^([0, +00)) with A strictly convex so that A'{p) is 
invcrtible. Further conditions on A{p) will be stipulated below. The first term in 
(1.1) represents then the internal energy of the star, the second term the rotational 
kinetic energy and the last term the gravitational potential energy. 

In [14] the existence of minimizcrs of J under the constraint (1.2) has been 
obtained. [13] contains further results for this model of prescribed angular velocity. 
In [11] support estimates for critical points of (1.1) under the constraint (1.2) have 
been shown. In particular, [11, Theorem 1] states that for u > coq > 0, the support 
of p is contained in a ball ^^(0, 0, ^3) for some ^3, where a = (7{u)o)- It follows that 

0<Bp<C in R^. 

Furthermore, [11, Theorem 2] shows that the number of connected components of 
the set {p > 0} is finite for any critical point p. 

Critical points of J with the constraint (1.2) are according to [11, (0.6)] charac- 
terized by the problem 

(1.3) p is continuous and nonncgative and A'{p) — ^uj^r^ — Bp = X{uj) in {p > 0}, 

where X{lli) is a Lagrange multiplier arising from the constraint (1.2). The focus in 
this paper is to study the free boundary d{p > 0} arising from (1.3). 

There is another model of rotating stars which has been studied in the literature, 
where the angular momentum per unit mass is prescribed. Existence of minimiz- 
ers for this alternative model has been obtained in [5], and the study of critical 
points has been carried out in [15]. Caffarelli-Friedman investigated in [7] the free 
boundary of minimizers for this alternative model. As Caffarelli-Friedman deal 
with minimizers, they are able to apply rearrangement methods to their functional 
to obtain solutions that are increasing in one direction which simplifies the analysis 
as well as the result. Unfortunately this technique does not work for critical points 
in either model and creates a difficulty for our analysis. Let us remark that the 
proofs presented in this paper for critical points of J with the constraint (1.2) work 
equally well for the study of the free boundary of critical points in the model in [7] . 
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The principal difficulty we encounter in our classification of singularities of the 
free boundary is that the nonlinearity is not an increasing function of the solution, 
so that various methods stemming from the well-known obstacle problem do not 
apply. Neither does the monotonicity formula derived in [1]. Let us also mention 
that our problem cannot be transformed into the type of problems studied in [8] , 
so we cannot use those results either. Another difficulty is that our equation is 
inhomogeneous. In particular, the leading order term on the right-hand side is not 
of the form f{u). This — together with a higher order degeneracy — distinguishes 
the present problem also from the recently researched "unstable obstacle problem" 
(see [16], [2], [4] and [3]). 

Last, let us point out that — due to the fact that the free boundary d{p > 
0} does not necessarily coincide completely with the free boundary of the PDE 
problem obtained by transformation — wc obtain in our classification of singularities 
several cases later called "pscudo cases." Wc suggest that in the case of minimizcrs, 
rearrangement techniques similar to those used in [7] may be used to show that 
solutions are decreasing in a certain direction, thus ruling out the pseudo cases. 

Setting u = \iJ^r'^ + Bp + A(a;), we obtain in the set {p > 0} that u = A'{p) = 
^~^{p), where $ : [0, -|-oo) ^ M is an increasing function satisfying according to 
the asymptotics 

A{p) = cipi + o{pi ) as p — )• 0, A{p) = C2p^ + o{pi) as p ^ +oo 

(where ci,C2 are positive constants) from Chandrasekhar's book [9, Chapter 10] 
and [11, (0.2)] the asymptotic relations 

(1.4) A'{p) = ^cipi -Fo(p3) as p -s- 0,A'(p) = ^C2ps +o(p3) as p ^ 0, 

(1.5) lim z-^/'^^(z) = c> 0. 
It follows (cf. [11, (3.3)]) that 

Au = 3w^ - 47rp = Sw^ - 47r$(u) in {p > 0} 

and 

Au = Sw^ in {p = 0}. 

Note that as u = A'{p) = ^~^{p) is only valid in the set {p > 0}, we obtain 
{p > 0} C {ti > 0} but not necessarily the opposite inclusion. It is however true 
that d{p > 0} = d{u > 0} and that if p(a;°) > then the connected component 
of {p > 0} containing coincides with the connected component of {u > 0} 
containing x^. 

Normalizing the equation as well as $ we obtain the free boundary problem 

(1.6) Au = 1 — $(u)xa with an open set ft satisfying xn < X{«>0}; 
the equation is to be understood in the sense of distributions. 




Figure 1 . Cusps (with the set painted red) 




Figure 2. Wedges (with the set Q painted red) 

Theorem A. For each solution u of problem (1.6) the following holds: Apart from 
the singular set S — {u ^ Qi\C] {Vu = 0} the level set {u = 0} and the boundary dfl 
are locally C'^'" -curves. The singular set contains in each bounded subset o/M^ at 
most finitely many singular points with the following possible asymptotics ( after 
rotation of the coordinate system): 

1. u{x^^ + tx)/t^ as t 0, and {u > 0} and fl satisfy the asymptotics 
demonstrated in Figure 1. 

2. There is X £ (-(X),0)U(1, +oo) such that u{x° +tx) / 1^ -> \ {Xxl + (1 - X)xl) 
as t 0. We may assume that X > 1 in which case {u > 0} and fl satisfy the 
asymptotics demonstrated in Figure 2. 
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3. There is X e (0, 1) such that u{x° + tx)/t^ \ (Axf + (1 - X)xl) as t 0. 
The complement of {u > 0} and that of fl is the single point x° . 

2. Proof of the Main Result 

Let u be a solution of (1.6). By LP- and C"-estimatcs u e W^^^{E.^) D Cl^"{R^) 
for each p S (1, +oo) and a G (0, 1). Differentiating u we obtain 

(2.1) AdkU = -^'iu)xndku € L^.iM.'') 

in the sense of distributions. Consequently u G W^iof (H^'^) ^ ^loc^^'^) ^"3^ each 
p € (1, +oo) and a G (0, 1). From (1.6) we infer now that the Hessian of u satisfies 

\D^u\ > c(n) > on {u = 0}. 

As {w = 0} n {Vm 0} is by the implicit function theorem locally a C^'"-surface 
— the regularity of the surface can be improved to real analyticity by the methods 
in [7] — , we will focus on the singular set S = {u = 0} Ci {Vm = 0}. From now on 
we will assume that u is axisymmetric, that is u = u{r,xs) and confine ourselves 
thus to a two-dimensional analysis. 

At each x^ G S we may rotate axes such that 

D'^u(x°) = I I for some A G ffi. 

^ ' \0 1-XJ 

Case 1: If < A < 1, then fi"^ consists in a sufficiently small ball Bs{x^) of only 
the point .t" which is in this case a local minimum point of u. 
Case 2: If A > 1 or A < 0, {m = 0} consists of two C^-curves intersecting at a 
nonzero angle at x^ (cf. Figure 3 and Figure 4): we may assume that x^ = and 
that A > 1. As in this case for sufficiently small 6, dx^u > in B^. n {xi > 3} and 
dx^u < in BrCi {xi < —S}, we may rescale and obtain that {u = 0} \ {0} consists 
of four C^'"-graphs. The fact that 

u{tx)/t'^ ^ {Xx\ + (1 - X)xl) ast^O 

implies now that the graphs have tangents as a; ^ and that we may combine 

them to two -curves intersecting at a nonzero angle at 0. 

Case 3: If A = 1 or A = 0, then {u = 0} consists either of two C^-curves ending 
in a cusp at x° (cf. Figure 6) or intersecting in a double cusp at x^ (cf. Figure 7): 
we may assume a;*^ = and A = 1, implying that u{tx)/t'^ \x\ as f — )• 0. Prom 
(2.1) with k = 2 we obtain that 

Ad2U = c{x)d2U 

with Holder continuous coefficients c{x). Applying [6, Lemma 3.1] repetitively for 
P = 3/2, 7/2, 11/2, 15/2, ... we infer that either 

d2U =p + T 
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n > 


\f u > 1 

>%\ 




n < ^ 



Figure 3. Double Wedge (with the set f2 painted red) 
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n < / 






/ 


n < ^ 



Figure 4. Pseudo Wedge (with the set Vl painted red) 

where p is a nontrivial harmonic polynomial of degree + 2 with leading term of 
order > 2 (by the fact that u{tx)/t^ — > as i — ^ 0) and 

|r(x)| < C^\x\P+^ 
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or d2U vanishes of infinite order at 0, that is 



\d2u{x)\ < Ck\x\'' in BrM 



for every k Cz N. In the latter case we obtain by repetitive apphcation of a weU- 
known strong unique continuation property (see [12, Remark 6.7] for a very general 
result), that 

d2U = in each connected component of {u ^ 0} touching the origin, 
implying by our information on the blow-up limit that 

u = + /(xi) in each connected component of {u ^ 0} touching the origin, 

where f{z) — > as z — > 0. But this contradicts the constraint J^s p = 1, thus 
proving that infinite order vanishing is not possible. 
Let us return to the former case 



where q is a nontrivial homogeneous harmonic polynomial of degree k > 2. It is 
important to note that by the fact that q is harmonic if (j((0, 1)) = then 
dxMiO, 1)) / 0- Similarly, if g((0, -1)) = then 5^,(z((0, -1)) ^ 0. 
It follows that 



In this case {u < 0} is in a neighborhood of a one-sided or two-sided cusp 
(depending on the signs of q{{0, 1)) and g((0, —1))); see Figure 6 and Figure 7. In 
the special case that q{{0,l)) = or q{{0,—l)) = 0, we obtain ^^^^^((0, 1)) ^ 
or dxiq{{0,—l)) ^ 0, respectively, in which case we obtain a non-symmetric cusp 
(even asymptotically) with respect to the xi-axis on the side of (0,1) or (0,-1), 
respectively. 
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Figure 6. Double Cusp (with the set ft painted red) 
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Figure 7. Cusp (with the set V, painted 



